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INCOMPRESSIBLE EULER EQUATIONS AND THE EFFECT OF 

CHANGES AT A DISTANCE 

ELAINE COZZI AND JAMES P. KELLIHER 


Abstract. Because pressure is determined globally for the incompressible Euler equa¬ 
tions, a localized change to the initial velocity will have an immediate effect throughout 
space. For solutions to be physically meaningful, one would expect such effects to de¬ 
crease with distance from the localized change, giving the solutions a type of stability. 
Indeed, this is the case for solutions having spatial decay, as can be easily shown. We 
consider the more difficult case of solutions lacking spatial decay, and show that such 
stability still holds, albeit in a somewhat weaker form. 


1. Introduction and Statement of the Main Theorem 
The Euler equations governing incompressible inviscid fluid flow in are given by 

{ d t u + u ■ Vw = — Vp 
div u = 0 
u\ t=0 = it 0 . 

There is extensive literature on existence and uniqueness of solutions to the Euler equa¬ 
tions on M. d in various function spaces. Short time existence of smooth solutions to (E d ) 
goes back to papers of Lichtenstein and Gunther [12, 9]; [10] contains a modern form of 
these results. For (E 2 ), surveys on global in time existence of solutions can be found in 
[13] and [4], A classical result for weak solutions of the Euler equations in the plane is 
due to Yudovich [20], which gives existence and uniqueness of solutions with velocity in 
L 2 (M 2 ) and vorticity (curl of the velocity) in L p D L°°(M. 2 ) for some p < oo. (Yudovich 
worked in a bounded domain, but his argument is easily extended to the full plane.) In 
recent years, some attention has also been paid to the Euler equations with initial vortic¬ 
ity which does not decay at infinity or may be unbounded (see, for example, [1], [7], [15], 
[16], [17], [18], [21]). 

We consider the question of what happens to the solution to the Euler Equations in W l 
at time t > 0 when a change is made to the initial velocity far away from a given point. 
We focus, in particular, on solutions to the Euler equations in two or three dimensions 
which have bounded vorticity and bounded velocity, but for which vorticity has no spatial 
decay. 

For solutions having spatial decay of vorticity, bounds on the effects of far-away changes 
can be obtained in an elementary manner in the L 2 -norm of the velocity. We quickly 
review the argument, as it helps motivate what follows. 
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Write B r (x) for the ball of radius R > 0 centered at x G R d . Let <p G Cc?(M d ) be 
supported in _B 2 (0) with ip = 1 on L>i(0), and dehne <Pr(x) = <p(x/R ). Let (ui,p±), 
(u 2 ,P 2 ) be two solutions to the Euler equations, and let w = u\ — U 2 , p = Pi — P 2 - Then 
subtracting (E d ) for u 2 from ( E d ) for ui, multiplying by (p R w, and integrating over time 
and space, one obtains 


( 1 . 1 ) 


IIw'WIIl* = I \<pr w ° 



2 (<p R Vip R ■ Ml) |w| 2 - 2 (ip R w • Vh 2 ) • (<p R w) 
+ 4 pwrVipr-w . 


For simplicity, suppose Vu 2 G L°°([0,f] x M d ). Then, assuming that our solutions have 
sufficient decay to ensure that 

2 


.4 11Vrt 2 1 1 ([ 0j *] X M d ) + 

3= 1 


M illL°°([0J]xR d ) T \\ U j\\L°°(0,t;L 2 ) " r IUb llL°°(0,t;L 2 ) 


+ I \Pj 


< oo, 


we have 


11 ^( 011 ^ < \\<Prw 


Oil 2 

II 2 


f*t " 


+ 


'0 L 


A 3 + 2A \\(p R w(s)\\ 2 L2 + 


4t1 2 

“r - 


ds. 


Gronwall’s lemma then gives the bound, 


( 1 . 2 ) 


\\ip R w(t)\\ 2 L2 < 


Wrw 


0 112 Cmax{l,4} 3 


lL 2 


+ 


R 


^2 At 


The bound in (1.2) controls how rapidly the values of the initial velocity outside of 
I?2 _r( 0) can affect the value of the velocity inside B R ( 0), showing that, for decaying solu¬ 
tions, the Euler equations have a degree of locality. (Of course, many variations on this 
theme are possible.) 

If we only have bounded, but decaying vorticity, say vorticity lying in L°°([0, t ]; L 1 flL 00 ), 
then we cannot conclude that Vu 2 is bounded. Regardless, an argument like that in [21], 
employing the Calderon-Zygmund inequality, can be used to handle the integral involving 
Vm 2 in (1.1). Osgood’s lemma is used in place of Gronwall’s lemma, leading to a different 
bound, but we suppress these details in this brief overview. 

Mere membership of tq, w 2 in L°°([0, £]; L 2 ) gives that \\(p R w(t)\\ L2 —* 0 as R —» oo; it 
does not, however, give a rate of convergence. The rate of convergence in (1.2) comes 
from bounding the rate of decay of the “tail” of the velocity, and through it the pressure. 

If we now assume that the velocity and vorticity are bounded but do not decay at 
infinity, the type of argument that led to (1.2) fails completely. The lack of decay in¬ 
troduces several difficulties when estimating terms in (1.1). First, while it is known that 
in 2D, Vp G L°°([0,f] x M 2 ) and \p(t,x)\ < C(t)log(|x| + e) ([15, 11]), this information 
is insufficient to control the pressure term in (1.1). Second, the integral involving u\ in 
(1.1) can no longer be controlled since w ^ L 2 (R d ). Third, the integral involving Vu 2 
cannot be controlled using the approach of Yudovich in [21], as the vorticity only lies 
in L°°(R a! ), and the Calderon-Zygmund theory does not apply in this setting. Finally, 
we would expect ||<Pa' u; (^)|Il 2 to grow with R anyway; hence, it would appear that such 
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an estimate is unattainable when the solutions do not decay. We should seek instead a 
weighted L°°-based estimate. 

One can view the bound in (1.2) as a type of stability estimate that demonstrates a 
partial locality of the solutions. Without some such locality, the very physical meaning of 
the solutions must be called into question: If a small change to the initial velocity at an 
arbitrarily large distance from a given point can immediately and significantly change the 
value of the solution at that point, then one could never effectively compute a solution. 

The original motivation for this paper was to address this issue of physical meaning 
for the case of bounded vorticity, bounded velocity (Serfati) solutions to the 2D Euler 
equations. Such solutions were first derived by Serfati in [15], where their existence 
and uniqueness was obtained. Serfati’s uniqueness argument easily extends to stability 
estimates in the L°°-norm of the velocity, but that argument cannot be localized. Locality 
of Serfati solutions has been an unaddressed problem. 

Stability estimates in L°°-based norms for Serfati solutions were established in [1, 17]. 
The estimate obtained in [1] is an outgrowth of Serfati’s uniqueness proof in [15], while 
the estimate in [17] is an outgrowth of Vishik’s uniqueness proof in [18], and both [15] 
and [18] depend critically on L°° estimates. (Neither employs an energy argument, nor 
is the maximum principle a key ingredient.) Our approach, which we describe in more 
detail below, will be to introduce a weighted L°°-norm into Vishik’s argument, adapting 
as well some pressure-based estimates in [17]. 

The weight we use is defined as follows. For fixed z EW 1 and a > 0, define g z : M. d —» M 
by 

9z(x) = (1 + \x- z\)~ a . 

A higher value of a makes the weighted norm more localized. As one might expect, the 
behavior of the Eulerian pressure limits how large a can be. In particular, the important 
properties of the Riesz operators described in Theorem 7 below play a crucial role in 
determining this upper limit. 

For a velocity held, u, on R d , we define the vorticity, co = oj(u), as 

u! = Vw — (Vu) T . 

(In 2D, one can alternately use u = diu 2 — d 2 u x .) 

We now state our main result. 

Theorem 1. ForfixedT > 0, assume that (u\,p\) and (u 2 ,p 2 ) satisfy (Ed) inD'((0,T) x 
M d ), with 

u u u 2 E L^dO^fLPnL 00 ^)) forp E [1, oo], 

Wl , W2 er([o,T] ; r(R rf )), 

d d 

Pi = RkRiu\u \, p 2 = RkRiu 2 u l 2 in BMO. 

k,l =1 k,l =1 

M = M(t ) := sup (1 + ||o;i(s)||l°° + 11^2(3)||l°° + ||rii(s)||£oo n £ P + ||tt 2 (s)||/ / ° o nLp) ? 

0<s<t 


Let 
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and set e 0 = ||skK - u° 2 ) || L°° • Define the strictly increasing function, p\ [0,1] ha [0,1] by 
p(x) = x(l — logx). Then, whenever a < 1 + d/p, 

\\g z (u! - u 2 )(t)\\ LOO < Mp (min | Cj^ + CMtp ((||) ) ,i}) • 

The constant C depends only upon d and a. 

Remark 2. We can see from Theorem 1 that for any t <T, we obtain a useful bound on 
\\g z {ui — u 2 )(t)\\ LOO as long ||^K-K)II Z/°° is sufficiently small. Also, stability in C' r (M d ) 
for all r < 1 follows from Theorem 1 by interpolation. 

We now briefly describe three examples to which Theorem 1 applies. All of these 
examples are 2D, for even though Theorem 1 is not dimension-dependent, only in 2D are 
weak solutions to the Euler equations known to exist. 

If p = oo and d = 2, then Theorem 1 above addresses Serfati (bounded vorticity, 
bounded velocity) solutions to the Euler equations in the plane (see [15]), and supplies 
an alternate proof of uniqueness to those in [15, 1, 17]. Moreover, it is shown in [8] that 
for Serfati solutions, 

II' u j(^)IIl°° — ciKlU°° (i + IK 9 IU~0 • 

Hence, M(t) < C(M( 0) + M(0) 2 t) for Serfati solutions. 

If u° belongs to L 2 (M 2 ) and cu° belongs to L°°(M 2 ), then there exists a global-in-time 
solution to (E 2 ) which conserves energy and for which vorticity remains bounded (see [6] 
for details). For such solutions, one can apply a simple argument using Littlewood-Paley 
theory to show that for all t > 0, 

\\u(t)\\ L °° < C(\\u(t )\\ L 2 + ||cn(i)||ioo) = C(\\u°\\ L2 + ||a; 0 ||i,oo). 

Thus, Theorem 1 applies for this class of solutions with p = 2 and d = 2, and we have 
M(t) < CM{ 0). 

If u° belongs to L 2 (R 2 ) and c o° belongs to L p D L°°(1R 2 ) for some p < 00 , then Theorem 
1 addresses the classical Yudovich solutions to ( E 2 ) and we again have M(t ) < CMfiS). 

Let us consider one simple application of Theorem 1, say, to be explicit, in the setting 
of Serfati solutions. Suppose that u\ = 011 -Br(O) for some R > 0 and that v!{ is 

compactly supported on B r( 0). Then 

(1.3) e 0 = \\do(Ui — U^Wl^ < ll' u 2lU° o (Bfl(0) c '(l + -R) °- 

The bound that results from Theorem 1 at time zero is 

|K - «2)(°)| < ll^oK - «S)|| £ . < CMp(Ce 0 /M). 

Because p is superlinear (linx r _ ) . 0 + p'(x) = +cxo), the bound at time zero is somewhat 
higher than the initial difference. Nonetheless, because the bound on eo in (1.3) decreases 
inversely with R , the bound on ||goK — u 2 )(t) || Loo decreases as R —* 00 for any fixed 
t > 0. This allows us to determine the distance, R , away from a given point (the origin, 
in this case) beyond which the initial velocity will not affect the value of the velocity at 
the given point up to a desired time. 

This example illustrates that the change to the initial velocity in Theorem 1 need not 
be localized, it need only be far away from the point or points of interest. We see this for 
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spatially decaying solutions as well in (1.2), though there the control is stronger and can 
be more cleanly expressed. 

Our proof of Theorem 1 is a localization of the uniqueness proof of Vishik’s in [18], 
with the pressure, however, handled much as in [17]. We localize Vishik’s argument using 
the function, g z . Each of the terms that appear in [18] have an analog here, though most 
of these terms must be bounded in a different manner. At the heart of bounding these 
localized terms are the inequality, 

(1-4) g z (x)h z (x-y) < (1 + \y\) a , 

established in the proof of Lemma 6, and the identity, 

(1-5) Vg z = i) z g zi 

where i\) z is a vector-valued function with L°°-norm independent of z G (in fact, 
lIVtzll l °° = «)• 

The identity in (1.4) is used several times, in each case to bound terms of the form 
g z (x)(<t> * f)(x), as follows: 

g z (x)(<f> * f){x) = j <j>{y)g s {x - y)g z {x)h z (x - y)f(x - y) dy 

< J lf%)l (1 + \y\r \(gj)(x - y)| dy = (|0(.)| (i + H)“»lfc/1)w. 

When (j) has sufficient decay, this effectively brings g z into the integral to obtain an 
estimate based on the function g z f. 

The identity in (1.5) is used to control the one term that is the analog of the integral 
involving u\ in (1.1). As a consequence of (1.5), \Vg z \ < C\g z \. Such a bound is never 
possible for a compactly supported function: this inability is the fundamental reason 
why we cannot use the “cleaner” localization that a compactly supported function would 
provide. This same issue shows up, for instance, in [8] (see Section 3.3), though it is 
addressed in a different manner. 

In [18], decay of vorticity is assumed, which allows the pressure to be easily dealt with. 
Without such decay, the pressure requires much more careful treatment. We follow the 
approach in [17], localizing it, much as we localize the approach in [18]. 

Remark 3. Being a localization of the approaches in [18, 17], we expect that a result 
similar to Theorem 1 will hold for vorticity in the spaces containing unbounded functions 
(6mo(R d ), for example) considered in [18, 17]. However, for simplicity, we do not consider 
unbounded vorticity here. 

Finally, we mention that the proof of uniqueness of Serfati solutions in [15] does not rely 
upon paradifferential calculus, and hence is not so tightly tied to the full plane. Indeed, 
Serfati’s approach was extended to a 2D exterior domain in [1], The proof involves 
bounding the difference between the flow maps for the two solutions (in this, it has 
similarities to the proof in [14] of uniqueness for Euler solutions in a 2D bounded domain). 
This approach seems resistant, however, to localization. It thus remains an interesting 
open problem whether an analog of Theorem 1 can be obtained for bounded vorticity, 
bounded velocity solutions in a 2D exterior domain. 
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2. Definitions and Preliminary Lemmas 


We first define the Littlewood-Paley operators. It is classical that there exists two 
functions x,</> G S(R d ) with supp y C {£ G R d : |£| < |} and supp </> C {£ G R d : | < 
|£| < |}, such that, if for every j > 0 we set fj(x) = 2 : > d (j)(2^x) then 

= x+Yl = L 

j> o j> o 

For n G Z dehne Xn G S^R^) in terms of its Fourier transform x n , where Xn satisfies 

Xn(0 = 

j<n 

for all £ G R fi . For / G ^'(R^) dehne the operator S n by 


Snf =Xn* f. 

Finally, for / G S'(R d ) and j > —1, dehne the Littlewood-Paley operators A j by 

'x * /» 3 — -1 

j> o. 

We will also need the paraproduct decomposition introduced by J.-M. Bony in [2]. We 
recall the definition of the paraproduct and remainder used in this decomposition. 


A jf = 


Definition. Dehne the paraproduct of two functions / and g by 

OO 

T f9 = A if A i9 = S i-if A j9- 

i,3 3=1 

i<3~ 2 

We use R(f, g ) to denote the remainder. R(f, g) is given by the following bilinear operator: 

R(f,g)= Y. A */A j9 . 

i,3 

Bony’s decomposition then gives 

f 9 = T f g + T g f + R(f,g). 

We will make use of Bernstein’s Lemma in what follows. A proof of the lemma can be 
found in [3], chapter 2. 

Lemma 4. (Bernstein’s Lemma) Let ry and ry satisfy 0 < ry < ly < oo, and let p and q 
satisfy 1 < p < q < oo. There exists a positive constant C such that for every integer k , 
if u belongs to L p (R d ), and supp u C -B(0,?yA) ; then 

(2.1) sup ||<9 “w||l*z < C k X k+d( 'p~^' > \\u\\L P . 

\a\=k 

Furthermore, if supp u C C(0, rA, r 2 A), then 

(2.2) C~ k X k \\u\\ LP < sup \\d a u\\ LP < C k \ k \\u\\ LP . 

\a\=k 

We will also make use of Osgood’s Lemma in the proof of Theorem 1. We refer the 
reader to [5] for a proof of the lemma. 
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Lemma 5. Let p be a measurable positive function, let 7 be a locally integrable positive 
function, and let p be a continuous increasing function. Assume that for some (3 > 0, the 
function p satisfies 

P(t) <P + [ 7(s)Mp(s)) ds. 

Jto 

Then -<j>(p(t)) + 0(0) < 7 (s) ds, where 0(x) = f* ^ dr. 

In what follows, we set 

h z (x) = 1/ g z (x) = (1 + \x - z\) a . 

The following lemma is a key technical ingredient in the proof of Theorem 1. 


Lemma 6. Assume f belongs to L°°(M. d ). There exists a constant C > 0, depending only 
on a, such that for every j > —1, 

(2-3) \\g z Ajf\\ L oo < C\\g z f\\ L ~, 

and, for 1 < k < d, 

(2.4) HsAAj/IU- < C2 j \\g z f\\ LO o. 

Proof. We first prove (2.3). Observe that for any a > 0, 

(2.5) *(«)*,(« -v)=( 0'7: y -J )"< ( 1+ . k r; l+ , l!<l )“<a+bi)« 


1 + \x — z 


1 + \x — z 


Assume j > 0. By (2.5), 

9z(x) / (j>j(y)f(x - y) dy = g z {x) / <t>j{y)h z {x - y)g z (x-y)f(x - y)dy 
jR d Jm d 

(2.6) < C\\g z f\\ L co f \(f>j(y)\g z (x)h z (x-y)dy<C\\g z f\\ L cof \(f>j(y)\(l + \y\) Q dy 

JRd J R d 

<C\\g,f\\ L - [ \(j)(y)\(l + \2~ j y\) a dy < C\\g z f\\ L °o. 

jR d 

An identical argument can be applied for j = — 1 with 0 replaced by y. This proves (2.3). 
To establish (2.4), note that for j > 0, 

g z (x)d k Ajf(x) = g z (x) / d k (p j (y)f(x - y) dy 

Jm d 

(2.7) =2 A d + 1 ) f (d k( j>)(2 J y)g z (x)h z {x-y)g z {x-y)f(x-y)dy 

JR d 

< C2 j \\g z f\\ L ~, [ (d k (/))(y)(l + \2~ j y\) a dy < C2 j \\g z f\\ L ~. 

J 

Again, an identical argument can be applied for j = — 1 with 0 replaced by X- This 
completes the proof of Lemma 6. □ 


In the proof of Theorem 1, the main restriction on the decay of g z , and hence on 
the value of a, comes from the pressure term. To estimate this term, we will need the 
following theorem from [19], which gives a pointwise bound on the Riesz operators applied 
to sufficiently decaying functions with moment conditions. 
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Theorem 7. Assume f is a function on which satisfies: 

1) \f(x)\<C(l + \x\)- d - N +% 

2) \D a f(x)\ < C( 1 + |x|) d N 1+e for all |a| = l, and 

3) J xPf(x) dx = 0 for all \(3\ < N 

for some fixed e G [0,1) and some positive integer N. Then for each i, 1 < i < d, 

\Rif(x)\<C(l + \x\)- d - N+e+s 
for every 5 satisfying 0 < 6 < 1 — e. 

We will apply Theorem 7 to both f(x) = Vy(x) and f(x) = D^fi(x), where x and 0 
are as above. In particular, we will make use of the following corollary to Theorem 7, 
which is proved in [7]. 

Lemma 8 . Assume g belongs to S(R d ). Given e > 0, there exists C > 0 such that for 
each i, j with 1 < i, j < d, 

\RiR j 'Vr)(x)\ < C{ 1 + |x|)- d - 1+e . 

Moreover, if belongs to 5(M d ) and has Fourier transform supported away from the 
origin, then given any 7 > 0, there exists C > 0 such that for each i, j with 1 < i, j < d, 

\RiRj1>(x)\ < (7(1 + |^|) -7 . 

We will also make frequent use of the following well-known lemma. A proof of the 
lemma can be found in [16]. 

Lemma 9. Assume u is a divergence-free vector field in LfifiW 1 ) with curl u = to. There 
exists C > 0 such that for every j > 0, 

(2.8) ||AjVti||£,oo A (7||Aju;|| L <x>. 

3. Proof of Theorem 1 


The proof of Theorem 1 is a modification of the proof of uniqueness of solutions to the 
Euler equations, as given in [18, 17]. The proof below deviates from those in [18, 17] in 
that we must account for the weight g z when proving the necessary estimates. To do this, 
we make use of Lemma 6 and Lemma 8 . 

Let «i, U 2 , g z , and T > 0 be as in the statement of Theorem 1. Set w = U\ — u 2 - For 
any N > 0, 


N 


\\gzw\\ L °° < II 9z&M\l°° + II 0 A 


w\\l° 


3=~ 1 


j>N 


< CN sup H^Ajinll 

L°° Eh Ajw\\ 

L°° 


-1 <j<N 


j>N 


< CN sup \\g z Ajw\\L°° 

j 

< CN sup \\g z Ajw\\L°° 

j 


+ c y 2 j (nAjOJiii^oo 

j>N 

+ C ^ 2 J (|| AjtuiUioo 

j>N 


+ 

+ 


II AjO^Hl 00 ) 
II AjU^IIloo) 


(3.1) 
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where we used Bernstein’s Lemma and Lemma 9 to get the third inequality. To estimate 
hzw\\ L°° 5 we establish an estimate for sup • \\g z Ajw\\L°° and only then judiciously choose 
a value of N (see (3.28)). 

We now estimate sup^ \\g z Ajiu\\L°°, a task that will consume most of our effort. 

Set p = pi — pi- Note that w satisfies 

d t w + w ■ Vw 2 + ui • Vw = —Vp. 

Applying the A j operator to this equation gives 

d t A jW + Sj_ 2 wi ■ VA jw = —Qj(ui, w) — Aj(w ■ Vn 2 ) — VA jp, 

where 

(3.2) Qj(ui,w ) = A j(ui ■ Vic) — S'j_ 2 ui • VA jw. 

We now multiply through by g z to get 

d t (g z Ajw) + g z Sj_ 2 ui ■ VA jw = -g z Qj{ui,w) - g z Aj(w ■ V« 2 ) - g z WApp. 
Observe that, by the product rule, 

g z Sj- 2 ui ■ VA jw = Sj_ 2 ui ■ gAJAjW 

= Sj- 2 ui ■ W(g z Ajw) - ( Sj- 2 ui ■ Vg z )AjW. 

This implies that 

(3.3) d t (g z Ajw) + Sj- 2 ui ■ V(g z Ajw) = / + // + III + IV. , 
where 

/ = {Sj-iUx ■ Wg z )Ajw 

II = -g z Qj(ui,w) 

III = -g z Aj(w ■ Vu 2 ) 

IV = -g z V A jp. 

Fix p 0 > 3 to be specihed later. We will estimate I-IV for three cases: 1) middle 
frequencies, 3 < j < p 0 ; 2) high frequencies, j > p 0 ; 3) low frequencies, — 1 < j < 2, in 
that order. 

So assume that 3 < j < p 0 . We begin with /. A simple calculation gives (1.5), which 
allows us to write 

II-^IIl 00 = \\(Sj_ 2 Ui ■ ij) z g z )Ajw\\Loo < a||5' j _ 2 Wi|U-||^A j u;|| i0 o 

(3 ‘ } < CM\\g z Ajw\\ L oo. 

Here, we used ||S' j _ 2 mi||l°° < C'||«i||loo. 

Next, we estimate III. Expanding III using Bony’s paraproduct decomposition gives 

III = -g z A 3 ^ Aiw ■ VSi_ 2 u 2 - g z Aj ^ S t _ 2 w ■ VA { u 2 
|i-j|<3 |K? I <3 

- 9zAj ^2 Aiw-VA^u-2. 

\l—i r | < 1 , 
max{l,i / }>j- 3 
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The first term in (3.5) can be estimated using Lemma 6 as follows: 

«A E Aiw ■ VS)_2'U2 E C \\g z Aiiv ■ X7 Si- 2 u 2 \\l°° 

( 3 . 5 ) Ml<3 \l-j |<3 

< c ^2 \\Si- 2 ^U 2 \\L^\\gzAiw\\ L oo < CM(j - 2)sup \\g z A q w\\ L ™. 

|i-j|<3 9 

To obtain the last inequality above, we applied Bernstein’s lemma and Lemma 9 to 
establish the series of estimates, 

1-2 

(3.7) ||S'z_ 2 Vm2||z,°° E || A_iVm 2 ||l°° + ^2 ||AjVw 2 ||.z,oo A C , ||m2||l°° + C(l — 2 )||cd2 ||z,°°- 

3=0 

For the second term in (3.5), another application of Lemma 6 and Lemma 9 gives 


g. 


■Aj ^2 Si- 2 w • VA|« 2 < C ^2 \\gzS 1 - 2 w ■ VAiM 2 || 


L°° 


(3.8) 


<c E (ii^vujIL- e II»a 

|Z—j|<3 V m<l—2 


M\l™ < CM{j - 2) sup \\g z A q w\\ L o 


The third term in (3.5) can be estimated using the divergence-free property of w, Lemma 
6 , Bernstein’s Lemma, and Lemma 9. We have the series of estimates 

d 

d k {Aiw k A l 'U 2 ) 

k= 1 


g z Aj ^2 A l wVAfU 2 

max{/,/ / }> t 7 ‘—3 

( 3 - 9 ) <C J2 2 J ||^A i w;A^ 2 || L oo < C 


| l-l | < 1 , 
max{/,/ / }>ji—3 


E 


2 j ~ l || A,/Vu 


l' V U 2 ||L C 


■Ilfl'zAziulIj 


\l — l | < 1 , 

max{Z,Z / }> 7 ‘—3 


| l-l |<i, 

max{Z,Z / }> € 7 ‘—3 


< C'||o; 2 || L <x> sup \\g z A q w\\ Loo < CM sup \\g z A q w\\ La °. 

Q Q 

Substituting (3.6), (3.8), and (3.9) into (3.5) yields the estimate 
(3.10) ||///||l~ < CM(j - 2) sup \\g z A q w\\ L oo. 
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We now estimate II. As in [18], we write Qj{u\,w ) = J2 m =iQT( Ul,w ^ w ^ere 


Q){ui,w) = Y, A jT dkW u}{ 

k =1 

d 

Qj(u 1 ,W) = -Yl T u* d k,A j]w 
k =1 
d 

Q 3 j(ui, w) = Y T v k 1 -S^ 2 uY A J W 


k =1 
d 


Q%Ui,w) = YA^j R ( u i d kw) - R(Sj- 2 Ui, Ajd k w)}. 


k=l 


To estimate g z Q](ui, w), we use the divergence-free property of u±, Lemma 6, Bernstein’s 
Lemma, and Lemma 9 to write 

d 


\g z Q)(u 1: w)\ = 


L°° 


g z Y A j®k Y Si-2U!^iUi Y ^WgzSi^wAiUiW, 

k=l |/-j|<3 |Z-j|<3 

< C Y - 2) sup \\g z A q w\\ L <x>\\A l Vui\\ L °a 

K-il<3 q 

< c 11 Will L oo (j - 2) sup \\g z A q w\\ L oo < CM{j - 2) sup \\g z A q w\\ L <*>. 

Q Q 

Next, we estimate g z Q 3 (ui,w). By properties of the Littlewood-Paley operators, there 
exists M 0 > 0 such that 

d oo 

Q$(u 1} w) = EE Si-2iu\ - S^uDAfikAjw 

k =1 1=1 


Y E ^- 2 E a ^ a ^ = E E E a ^a 


w. 


k =1 |i—j|<3 \m=j —1 


k =1 |i—j|<3 \\m—j\<Mo 


Hence, 


(3.11) 


7//, Vi.'.' 1 .) < Y E I E 2 m ||A m Vui||£,o. ) ||g,9 l A 1 A,TO|| I o 

k=1 \l~j l<3 \\m-j\<M 0 


< cilcuiiiioo y ^ y ~] 2 ^ 11 ii^00, 

fc=l \l-j |<3 

where Bernstein’s Lemma was used to get the first inequality and Lemma 9 was applied 
to get the second inequality. Now note that by Lemma 6, 
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Inserting this estimate into (3.11) gives 

\g z Q i j (u 1 ,w)\ <CMswp\\g z A q w\\ L oo. 

q 

To estimate g z Q?(ui,w), we observe that there exists M 0 > 0 such that 

d 

Qtju^w) = - y y [Si^u^A^dkAtw. 

k=l \l-j\<M 0 


We express the commutator as a convolution integral and utilize the divergence-free prop¬ 
erty of Mi to write, 

d 


Qj(u 1 ,w)(x) = — y y/ 4>j{y){Si- 2 u l i(x) - Si_ 2 Ui{x - y^dkAiW^x - y)dy 
\l-j\<M 0 k =1 ^ Rd 

d n 

= - y y / d k (j) j (y)(Si_2Ui(x) - Si_ 2 Ui(x -y))Aiw(x -y)dy. 

-iKMn fc=l ' dRd 


|J-j|<Afo fc=l 

Hence, 

< 0 *(a:) y 2 : j(d+1) / |V^(2- 7 i/)|||5' i _2VM 1 || L oo|i/|^(x-i/)/i 2 (x-i/)|A Z M;(a;-i/)|(ii/ 

|Z-J|<M 0 ^ Rd 

< y 2 j(d+1) ||5f z A z w|| L oo||S' z _2Vn 1 || i oc f \V(f){ 2 3 y)\\y\g z (x)h z (x-y)dy 

\l-j\<M 0 ^ Rd 

< y 2^'||^A Z M;|| i0 c||5 z _2VM 1 || L o 0 [ |V0(r/)||2-^|(l+|2-^|) Q di/ 

|i-i|<Afo ^ Rd 

< CM(j - 2) sup 11 sc Agio 11 z/x,, 

q 


where we applied (2.5) and a change of variables to get the second-to-last inequality and 
(3.7) to get the last inequality above. 

For g z Qj(ui,w), we again follow [18] and write Qj(ui,w) as the sum of two terms: 

Qjiu^w) = y i ('u i , w) + Q^iu^w), 


d 

Q^iu^w) = y A jd k R(u\ - Sj- 2 Ui,w), and 

k =1 
d 

Qf 2 (ui,w) = y {A j i?(6 , J _ 2 Mi, clfcw) - R(Sj- 2 Ui, Ajd k w)}. 

k =1 


where 
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By properties of Littlewood Paley operators and Lemma 6, 

d 

\g z Qf\u ll w)\= g z ^2Ajd k E A i( u i ~ S j-^i)A l 


>W 


k =1 


| i-i |<i, 
max{/,/ / }> < 7 — 3 


<£?E E yWgzMul-Sj^uDA^wWL- 


k — 1 |l—l' |<i, 

max{Z,Z / }> < 7 '—3 

< Csup \\g z A q w\\ L oo E 2 J ~ l \\A l Vu 1 \\ L ™ < CMsup \\g z A q w\\ L ™, 

q l>j-Mo q 

where Lemma 9 was again used to get the last inequality. 

The strategy for estimating Q*’ 2 (ui,w) is similar to that used to bound g z Q 2 (ui,w). 
By properties of Littlewood Paley operators and the divergence-free property of rq, 

d j —i ;+i 

Qf 2 {ui,w) = c^E E E [ A j,^iS j -2Ui]A l >d k w 

k=l l=j-M 0 

d j 1 Z+l t 

j- 2 u\{x) - AiSj^ 2 u\{x - y)\A l >d k w(x - y) dy 


EE E 2ii [ 

k=n= j -Moi:= l -i jRd 
d j —1 l+l „ 

EE E 2iW1) / (+m , v)\+s 1 

k=ll=j~Moi'=i-i jRd 


- 2 u\{x) - AiSj- 2 Ui(x - y)]A l 'w(x - y) dy. 


Hence, 

\g z Qf 2 {u 1 ,w)\ 


3-1 l+l 


< 


E E 2j(d+1) / \V(j)(2 J y)\\\Ai\/S j _ 2 Ui\\L^\y\gz{x)h z (x ~ y^lg^wWLoo dy 
“L,/, J* d 


l=j-M 0 i' =i—i 
j —1 Z+l p 

* E E 2^ +1 )||A i Vn 1 || i oc||^A^II^ / |V0(2L/)|M(l-H2/irch/ 

Z=j-Mo ■ /Rd 

i z+i /» 

^E E ||cji|| LC x>||^A;/w;|| l ^ / |V0(2/)|M(1+|2 J 2/|)"ch/ 

«=i-M 0 e=z~i jRd 

< CM sup \\g z A q w\\L°°, 
g 

where we again used (2.5) and Lemma 9. 

We conclude that 

\\II\\ L «> < CM(j - 2) sup \\g z A q w\\ L ™. 
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We now estimate IV. We follow arguments from [17]. Set (ftj = 4>j-\ + <f>j + <f>j+ 1, so 
that A j = Aj_i + A j + Aj + i. As in [17], write 

d 

IV = —g z VAjp = -^V(-A)" 1 ^ dkAfiAjiw'ul) 

(3.12) 

-^V(-A)" 1 ^ d^d.A j (w k u l 2 )=:V + VI. 

k,l=—1 

For V, observe that 

d 

V = g z {x)^2 / V{-A)- 1 d k <t> j (y)diA j (w l u'[)(x - y)dy, 


k,l=1 ' 


so that 


(3.13) 


ID < V ||s«5iA,'(to'ti{)|| L » / |V(—A) l dijj(y)\g,(x)hi(x-y)dy 

k,l= 1 


< 


II^Aj-^ii^Hioo / |v(—A) 1 9 fc 0 i (2/)l(i + M)“dy, 


fc,«=i 


where we used (2.5) to get the last inequality. The operator V(—A)~ 1 d/ C can be expressed 
as a product of Riesz operators, which commute with dilations. Thus, a change of variables 
gives 

d 

ID < J2 / |V(-A)- I a fc ^(y)|(l + |2-’3/ir d*/ 


(3.14) 


k,l=1 ^ K 

d „ 

< \\9zdiAj{w l u\)\\ L ^ / V(-A)” 1 40(|/)|(l + 12/1)“% 

k,l=1 l ' Rd 


Observe that, since 0 has Fourier transform supported away from the origin, Lemma 8 
can be applied. Thus the integral in (3.14) is finite, and 

d 

(3.15) \V\<CY t \\9zdiA j {w l u k 1 )\\ Loo . 

k,l=1 

We estimate VI in exactly the same way as V, and we conclude that 

d 

(3.16) \VI\ |M fc A> fe 4)|| ioo . 

k,l= 1 

We can now complete the estimate for V and VI by applying an argument identical to 
that used to estimate III. This gives 

II^IIlcc < ||R|| L cc + \\VI\\ LOO < CM(j - 2) sup \\g z A q w\\ L ~. 
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Viewing (3.3) as a transport equation (g z Ajw transported by Sj_ 2 Wi), integrating in 
time, and applying the estimates for I, II, III, and IV above, we conclude that for 
3 < j < Po, 

\\gzAjw(t)\\ L °o < \\g z A j w°\\ L ™ + CM -2)sup||^A q w;(s)|| L oo^ ds, 

so that 

(3.17) sup \\g z Ajw(t)\\ L oo < sup \\g z A j w°\\ L °o+CM / ( Po sup ||^ 2 A g in(s)|| I ,<x= ) ds. 


3<j<po 


3<j<Po 


For the case j > po, note that by Bernstein’s lemma and Lemma 9, 


(3.18) 


sup \\g z Ajw(t)\\ L °° < sup 2 J ||AjVw(t)||ioo < CM2 po . 

j>PO j>PO 


Now assume —1 < j < 2. In this case, ||/||l°° can be handled exactly as in (3.4). For 
III, the divergence-free assumption on iq and Lemma 6 imply that 

\III\ = \g z AjV ■ (wu 2 ) | < C2 j \\g z wu2\\L°° < CM\\g z w\\ L °°- 


For IV, we write 
(3.19) 


; VA jp - ' 7 . W RiRkVAj(w l u* + u\ 


w 


i,k =1 


We focus on g z R l R,k'VA 3 ( w l u \). Note that for 0 < j < 2, 
g z R i R k X7A j (w l u l l) = / (R i R k \ / (p j (y))g z (x)h z (x - y)v![{x - y)(g z w l )(x - y) dy, 


so 


(3.20) 


\g z RiR k WA j (w t u l l)\ < C'||ui||l<x.||^k;|| L oo / \RiR k V<j)j(y)\(l + \y\) a dy 

J R d 

<CM\\g z w\\ L oo [ \RiR k V(j)(y)\(l + \2~ j y\) a dy, 


where we used (2.5) to get the first inequality. Since V0 has Fourier support away from 
the origin, Lemma 8 can be applied to conclude that the integral in (3.20) is finite. Thus, 

\g z RiRCVA j {w l v![)\ < CM\\g z w\\ L ™. 


For j = —1, we must estimate 


(3.21) 


= g,(x) I R i R k X'x(y)(w'u1)(x - y) dy. 


i„,k\ 


This term is the only term which places a substantial restriction on the growth of h z (or 
decay of g z ) at infinity. Indeed, estimates for all other terms in the proof of Theorem 1 
will hold with h z {x) = (1 + \x — z|) Q for any a > 0. It is in the estimate for (3.21) where 
we utilize the membership of u k and u 2 in L p (R d ) to optimize the value of a. 

Let q satisfy l/p+ 1 /q = 1. Then 

g z R i R k \7A_i(w l u’l) = g z (x) / RiR k Vx(y)h z (x - y)u\{x - y){g z w l ){x - y) dy, 

J]SL d 
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and hence, 

\g z RiRk ^< || RiR k \7x(-)gz(x)h z {x - ■)\\ L q\\ui(g z w' l )\\ LP 

< \\R i R k Vx(-)(l + I ' 

where we used (2.5) to get the last inequality. 

Given p , we must determine values of a for which ||i?ji4,Vy(-)(l + | • |)“||l 9 is finite. In 
view of Lemma 8 , a short calculation shows that ||i?ji?fcVy(-)(l + | • |)“||l 9 is finite if 

a < 1 + d/p. 

Therefore, assuming a < 1 + d/p, it follows that 

\g z RiR k V< CM\\g z w\\ L oo. 

Of course, arguments identical to those above imply that 

\g z RiR h ^^i(u l 2 w k )\ <CM\\g z w\\ L ™, and 
\g z RiR k \/Aj(u l 2 w k )\ < CM\\g z w\\ L oo 

for 0 < j < 2. Inserting these estimates into (3.19) gives the following estimate for 
-1 <j< 2 : 

\g z X7Ajp\ < CM\\g z w\\ L °o. 

We now turn our attention to II. Instead of considering the commutator, we estimate 
each term in the difference individually. For the first term in the difference, we observe 
that, by Lemma 6 , 

\g z Sj- 2 Ui • VAjiu| < ||S' i _ 2 ni5' z VA j u;|| L oo 

< \\Sj- 2 U 1 \\ L °°\\gz'VAjw\\ L °o < CM2 J \\g z Ajw\\ L oo 

for — 1 < j < 2. For the second term in the difference, note that by the divergence-free 
property of U\ and Lemma 6 , 

\g z Aj(ui ■ Vw)| < C2 J \\u 1 g z w\\ L °° < CM\\g z w\\ L °o. 

We also recall the following estimate, which is derived in (3.1) and which holds for any 
Po > 0: 

(3.23) H^wUloo < Gposup \\g z A q w\\ L °o + CM2~ P0 . 

g 

Integrating (3.3) in time, applying the above estimates for I, II, III, and IV when 
—1 < j < 2, applying (3.23), and taking the supremum over — 1 < j <2 gives 

sup || 0 ~Ajw(t)||z,oo < sup II^Aj^Hloo + CM / (p 0 sup \\g z A q w(s)\\ L °° + 2~ po ) ds. 
-l<i<2 -1<J<2 Jo q 
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Combining this estimate with (3.17) and (3.18) yields 



sup||^Agiu(t)|| L <*, 

<7 



(3.24) 

< sup \\g z A q w°\\ 
<? 

Lea + C M j | 

Po sup \\g z A q w(s)\\ L ™ 
v q 


< Ce 0 + CM [ 

tf / 

(Vo sup \\g z A q 

rw(s)|U°° + 2 -P0 ) ds , 


Jo 

\ <? 

J 

where 






eo := \\gzW°\ 

1 L°° ’ 


>—p 0 


ds 


and where in the last inequality we applied Lemma 6. 

At this point, we could apply Gronwall’s estimate to obtain 

sup \\g z A q w(t)\\ Lao < (Ce 0 + CM2~ po t ) e CMpot 
q 

and choose p 0 in a manner so as to minimize the right-hand side. This will not, however, 
yield a sufficiently tight estimate. Instead, we will choose po in (3.24) to vary with t (this 
is valid, because the constant, C, in (3.24) does not depend on p 0 ) in a way that allows 
us to apply Osgood’s Lemma. 

Toward this end, first observe that it follows from Lemma 6 that for all t > 0, 


m := 


sup \\g z A k w(s)\\Lc 0 ds < 1, 


C^Adt J o k 
where Cq is the constant in (2.3). Now choose 

Po = 3 - log 2 S(t). 

Note that po > 3 as required. 

Substituting our chosen value for p 0 into (3.24) yields 

sup \\g z A q w(t)\\L°° < Coeo + CMt5(t ) + CM 2 t{ 3 — log 2 6(t))6(t) 

(3.25) 9 

< C 0 e 0 + CM 2 t5(t)( 3 — log 2 S(t)) < C 0 e 0 + CM 2 t5(t)( 3 — logh(f)). 


Setting t — s in (3.25) and integrating both sides from 0 to t, then dividing both sides 
by C 0 Mt gives 

m < + CMj / ®5(s)( 3 - logiW) ds < A: + CM fs(s)( 3 - log 6(s))ds. 

Since x ha x( 3 — logx) is increasing on (0,1] and we know that S(s ) G (0,1], we can 
apply Osgood’s Lemma to obtain 

- log(3 - log<S(t)) + log ^3 - log (-^)) < CMt. 

Taking the exponential of both sides twice, we conclude that 
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Now increase the bound in (3.25) to 

(3.26) sup \\g z A q w(t)\\ L oo < C 0 e 0 + CM 2 t5(t )(4 - log 5(t)) 

q 

and note that 5(t) < e 3 . Then since x ha x( 4 — logo;) is increasing for x G (0, e 3 ], we can 
replace 5(t) in (3.26) with S(t). This gives, 

sup \\g z A q w(t)\\ L oo < C 0 e 0 + CM 2 te 3 ^4 - log ^e 3 


(3.27) 

By (3.1), 


= C 0 e 0 + CM 2 te 3 (|| 


1 — e 


- C “> og (^)):=a 


■>-N\ 


11 11 — C(Na + MT 

Define F(N) := Na + M2~ N . This has a minimum when 
(3.28) IV = IVo := log 2 (Mlog2/a), 

where we treat N 0 as though it were an integer. At this (near) minimum, assuming that 

I L° 

Ma a a , a 


0 < a < M log 2, which holds for all sufficiently small ||(yf,w°|| f/X) , we have 


F(N 0 ) = a log 2 + 


= M- 


a M log 2 M log 2 
= M(b — (log2)61og 2 b) = Mb{ 1 — log b), 


— M log 2-log 9 - 

M log 2 2 M log 2 


where b := aj{M log 2) < 1. (Note that if a > M log 2, we can do no better than setting 
N = 0, which gives the useless bound, ||g 2 w|| LOO < CM.) 

Hence, when a < M log 2, we have 

\\g z w\\ Loa < CMb( 1 - log b), 

which is the bound stated in Theorem 1. 
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